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Abstract. A self-consistent fluid model developed for simulations of micro- 
gravity dusty plasma experiments has for the first time been used to model 
asymmetric dusty plasma experiments in a modified GEC reference cell with 
gravity. The numerical results are directly compared with experimental data 
and the experimentally determined dependence of global discharge parameters 
on the applied driving potential and neutral gas pressure is found to be well 
matched by the model. The local profiles important for dust particle transport 
are studied and compared with experimentally determined profiles. The radial 
forces in the midplane are presented for the different discharge settings. The 
differences between the results obtained in the modified GEC cell and the results 
first reported for the original GEC reference cell are pointed out. 
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1. Motivation for this study 

The GEC reference cell was originally designed to allow fair comparison between 
plasma processing studies performed in different laboratories [HE]. A large 
experimental and numerical effort was undertaken to understand both the proper 
operation of the cell from a technical viewpoint, as well as the physics of chemically 
active plasmas with different gas mixtures [3j |H [5] . 

In the early nineties, it was realized that the GEC cell could also play the 
role of a standard experimental platform for dusty plasma experiments [B]. Several 
modifications to the original design were required in order to suspend dust particles in 
the discharge and to use different optical systems to visualize them. This in a way led 
to a loss of standardization, since different solutions were to deal with the additional 
challenge of conducting dusty plasma experiments [7] . 

Despite numerous experimental and numerical efforts to describe these 
experiments, a study employing a self-consistent numerical model, which can self- 
consistently calculate the dust forces from the plasma parameters, is missing. 
Furthermore, the changes in discharge characteristics with respect to the original 
GEC reference cell, due to the necessary modifications for dusty plasma experiments, 
are usually ignored. 

We have developed such a code in the past and applied it to micro-gravity dusty 
plasma experiments in symmetrically driven discharges as well as to devices including 
the effects of gravity and additional thermophoretic forces, due to heated surfaces 
[H [9l [10l [11] . The results have always shown excellent agreement with results reported 
in the literature, but the model has never been directly compared to measurements in 
a GEC reference cell. 

The motivation for this study is therefore to examine dusty plasma environments 
in a modified GEC cell with a self-consistent dusty plasma model for the first time, 
to compare results from the model to measurements of plasma properties in the 
experiment, rather than from the literature alone, and determine the effect of the 
modifications to the GEC cell on the local and global discharge characteristics. The 
latter depend on the global particle and power balance of the discharge, which can 
be observed through the DC bias on the powered electrode and the power absorbed 
in the plasma. The local parameters include the dust charge, the plasma densities, 
and the plasma potential, which directly determine the forces that would act on dust 
particles present in the discharge. 

Section [2] describes the geometry of the modified GEC cell used at CASPER, and 
clarifies the Langmuir probe measurements, section [3] discusses the numerical model, 
section 2] shows the results for the global parameters studied, i.e. the DC bias and 
the absorbed power and section [5] shows the results for the local plasma profiles, i.e. 
the plasma potential, the particle densities, and the derived dust charge number. In 
section [6j the forces obtained from the measured plasma parameters are presented and 
compared to the outcome of the model. Section [7] discusses the results and briefly 
mentions the outlook for future work. 
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2. Description of the experimental setup 

2.1. The modified GEC cell. 

The GEC reference cell used by CASPER, shown schematically in figure [TJ is modified 
to allow dusty plasma experiments to be performed. The upper electrode is a grounded 
hollow cylinder rather than a solid electrode, so that a top-mounted camera can be 
used to take pictures of dust crystals from above. In total, two camera/dye laser 
systems have been added, to capture side-view and top-view pictures of dust clouds 
suspended in the plasma volume. The lasers are equipped with cylindrical lenses, 
to creat thin laser sheets that illuminate selected areas within the dust clouds. The 
cameras can also be equipped with filters that allow only light at the laser frequency 
to pass through. This helps to select light scattered by the dust particles, and not 
from the plasma glow. 

The dust clouds suspended in the modified GEC cell are confined in the radial 
direction by a parabolically shaped electric potential, created by a circular cutout 
milled in a cover plate set on top of the powered electrode. The different cutouts used 
in the experiments have radii of 0.63, 1.25, and 2.5 cm. To introduce particles into the 
plasma, two dust shakers have been added to the top flange near the upper grounded 
electrode. Tapping these shakers forces dust particles to enter the plasma under the 
force of gravity. The bottom of each dust shakers is covered by a calibrated mesh to 
prevent larger clumps from entering the plasma. 




Setup of CASPER's GEC RF Reference Cell. 

Figure 1. Sketch of the interior of the modified GEC cell. The lower electrode 
is powered, the upper electrode is grounded, as are the outer walls and the 
groundshield surrounding the lower electrode. The inter-electrode spacing is 3 
cm, the radius of the cell is 13 cm, the height is 21 cm, and the electrode radius is 
5.4 cm. The cutouts used for this experiment were 0.63, 1.25, or 2.5 cm in radius 
and 1 mm deep. 

The discharge parameters under external control include the neutral gas pressure, 
which can be adjusted using a butterfly valve controlling the input gas flow, the input 
power, which is adjusted by changing the driving potential of the RF source, and the 
DC bias on the powered electrode. For the current study, the DC bias was allowed to 
float. 
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2.2. Langmuir probe. 

The plasma parameters in the modified GEC cell described above were measured with 
a SmartProbe produced by Scientific Instruments LTD [12]. The SmartProbe is a 
RF-compensated Langmuir probe, with a 10 mm long, 0.38 mm diameter tungsten 
probe tip attached to a 470 mm long shaft. By moving the shaft toward or away from 
the center of the powered lower electrode, radial profiles of the plasma parameters 
were obtained at the midplane between the two electrodes. The probe was inserted 
into the system through a sideport in an attempt to ensure that the profiles were 
obtained exactly in this plane; however, bending of the probeshaft by gravity can not 
be excluded. The deviation from this plane was estimated to be less than 10% of the 
inter-electrode spacing, after inspection of side-view pictures. 

The radial probe position with respect to the electrode center was determined 
through examination of still-frame images from the top-mounted camera. Using 
imaging software [13j . the midpoint of the 10 mm long probe tip and the center 
of the electrode were determined. The line connecting these two points was taken 
to be the true radial direction and any small angle the tip made with this line was 
measured. The projection of the probe tip onto this line was then taken as the error 
in the radial direction, assuming that the probe measurements represent the plasma 
parameters averaged over the length of the tip. 

For each radial position, the bias voltage on the probe tip was swept from -95 V 
to +95 V in steps of 0.1 V. For every voltage step, the current collected by the probe 
tip was measured ten times, and the average was then computed by the probe software 
and stored. This measurement was repeated several times, for each radial position. 
Employing the IV-characteristic data, the probe software computed several plasma 
parameters using standard Laframboise theory for cylindrical probes [14] . Parameters 
computed include the electron and ion density, n e , rt+, the plasma potential, Vp, the 
floating potential Vfi, and the electron temperature T e . 

3. Description of the numerical model 

A two-dimensional hydrodynamic model is used to solve the equations for the electron- 
and ion-fluid (in an argon discharge), coupled to the dust fluid. In this study, we are 
concerned with the effect of changing discharge settings on the dust transport, and the 
effect of the modifications in the modified GEC cell on the plasma parameters. We 
therefore only consider situations where small numbers of dust particles are present in 
the discharge, so that the plasma parameters, and the forces acting on the dust, are 
not altered by the dust itself. However, the forces are still self-consistenly computed 
from the plasma parameters. We now proceed with a description of the solution for 
the plasma parameters. 

3.1. Solution scheme for the plasma. 

The continuity equation for the density nj for the electrons and ions (j = e, ions) is 
solved using a drift-diffusion approximation for the flux, Tji 

dn ' 

^ + v-r J = 5„ r »,i>,i<: n.x,,,. (1) 

with fij the mobility and Dj the diffusion coefficient. The sinks and sources Sj include 
electron-impact ionization and electron-impact excitation. The electric field is found 
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from the Poisson equation, 



v 2 y 



- (n+ - rie) , E = -VV, 



(2) 



with n+ , n e the ion and electron density, eo the permittivity of vacuum, and e the 
electron charge. 

Since the argon ions are too massive to follow the instantaneous electric field, 
an effective electric field is calculated to include the effect of ion inertia by solving 
dE e ff/dt = v m _ + (E — E e ff), with v m ,+ the momentum transfer frequency for ion- 
neutral collisions. 

A similar set of equations is solved for the average electron energy density, 
w = n e e, with e the average electron energy. 



In the above, — eT e ■ E is the Ohmic electron heating, which is the power input. The 
sinks, S w , include electron impact ionization and excitation. The ions are assumed to 
locally dissipate their energy, so that it is not necessary to solve a similar equation for 
the ions. Equations (jTJ [2J [3j) are progressed in time on sub-RF timescales until the 
solution set U(t) becomes periodic over a RF cycle to within a very small user-defined 
parameter; (U(t) = U(t + trf))- 

3.2. Solution scheme for the dust fluid. 

3.2.1. Dust particle charging A spherical dust particle with radius a immersed in 
plasma absorbs electrons and ions (with mass m e and m + respectively) until in 
equilibrium the electron and ion currents balance. Due to the high electron mobility 
compared to ion mobility, the equilibrium dust charge becomes negative, V(a) < 0. 
Using energy and angular momentum conservation, the ion and electron current can 
be calculated from OML theory [15] through 



For typical dusty plasma experiments, charge-exchange collisions increase the ion 
current to the dust particles, making the equilibrium dust charge less negative [16] . 
This is included in the above equations as the final term between square brackets in the 
ion current equation. In the equations, T e is the electron temperature, k is Boltzmann's 
constant, Xd is the linearized Debye length, and l m fp is the ion- neutral collision mean 
free path. E + = ^m + v 2 = 4 T+ + ^m+u 2 ^ is the mean ion energy, with T + ps T gas 
the ion temperature [17] and w+ is the ion velocity found from the flux, equation |T]) . 
The dust charge number is found using a capacitor model: eZo = ^toaV{a). 



— + V ■ T w = -eT e -E + S, 



W 1 
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3.2.2. Dust particle transport The forces acting on dust particles include gravity, the 
electrostatic force, ion drag, neutral drag, and thermophoresis. The gravitional and 
electrostatic force are calculated from 

F g = m D g = — —a 3 pg e z , F E = eZ D E, (6) 

where p = 1510 kg m~ 3 for melamine-formaldehyde particles, g = 9.81 m s -2 , and 
E the time-averaged electric field. The neutral drag and thermophoretic force can be 
calcuated using 

„ 4tt 2 32 2 K T VT gas 

F n = —a PgasVgasVD, F th = -— O , (7) 

G 10 Vg as 

with Vgas defined as the thermal velocity of the neutral gas, p gas — mA r n ga s the neutral 
gas mass density for argon, VT gas the temperature gradient of the background gas, 
and kt = 0.0177 W K _1 m _1 the thermal conductivity of argon at 300 K. 

The ion drag force is calculated including the effect of ion scattering beyond the 
Debye length, the effect of charge-exchange collisions and the effect of ion drift through 



Fid = n + m + vu + I a c (v) + np c {v) 



A 



D 



mfp 



(8) 



where v is determined from the mean ion energy, a c (v) is the OML collection cross 
section [TS], pc{v) the Coulomb radius, A(u) the Coulomb logarithm for scattering 
beyond the Debye length [18] . and v is found from a fit for the total energy in the 
Coulomb logarithm in order to take the effect of significant ion drift on the screening 
into account [TS]. K{x) = x arctan(x) + (•y/f — l) ~ \f\ m (l + x2 ) is the collision 
operator, used to calculate the effect of charge-exchange collisions on the ion drag |20j . 

Assuming a force balance between the neutral drag and all other forces, a drift- 
diffusion type equation can be derived for the dust flux, 

T D = nuimo^nDy 1 [F g + F E + F th + F ld ] . (9) 

with v m D the dust-neutral momentum transfer frequency. This equation can be used 
to transport the dust in time. In this study, we are simply concerned with the 
differences in the net force, due to changes in the geometry of the modified GEC 
cell, as well as the effect of the discharge settings on the forces. 



3.2.3. Geometry The cylindrical symmetry allows us to model half of the volume, 
approximated by a half-iJ-shape, and the two spatial directions become the radial 
direction, with coordinate r and the axial direction, with coordinate z. The biggest 
difference between the real and the model geometry is that in the model the top 
(grounded) electrode is solid, whereas in the experiment it is a hollow cylinder. The 
radial cutout plate, which provides radial confinement for the dust particles, is added 
in the model as an additional boundary condition for the potential on the inner part of 
the lower electrode, varying quadratically from the center of the electrode to the edge 
of the cutout. The numerical grid solution is 48 radial by 96 axial points. The vertical 
grid points are divided in three regions. Region I between the electrodes (2/3 of the 
points, giving a resolution of 3 cm/64 points = 0.5 mm/grid interval) and regions II 
and III next to the electrodes (both with 1/6 of the grid points). 

Once the plasma parameters are converged, the solutions on the grid are used to 
compute the dust charge, using the OML charging currents from equations |4] and 5. 
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Once the dust charge is determined, the electrostatic force and the ion drag force are 
computed. The neutral gas temperature is solved with an energy conservation scheme, 
which in this study only involves the heating by local energy deposition of ions and 
the thermal conduction by the gas. Once the temperature gradient is resolved, the 
thermophoretic force can be calculated. The sum of these forces thus gives the net 
force that would act on dust particles suspended in the plasma. 

4. Results for global discharge parameters 

The global characteristics can be measured by the absorbed power and the natural 
DC bias on the powered electrode. If our model is to correctly represent the actual 
experiment, the DC bias and the total power absorbed in the discharge should be 
similar to the experiment for given discharge settings. We therefore first compare 
these two easily determined global parameters. 

4--1- Natural DC bias. 

In asymmetric discharges, such as those found in the GEC cell, the number of electrons 
collected by the powered electrode per RF cycle during the positive phase of the cycle 
is unequal to the amount of ions collected during the negative phase. Therefore, 
the surface of the powered electrode acquires a negative charge, corresponding to a 
negative potential, called the (natural) DC bias. As a first approximation [55], it 
only depends on the relative surface of the plasma facing grounded area (A g ) and the 
plasma facing powered area {A^), and the amplitude of the driving potential, or the 
root-mean-square value, V rms — Vq/V2, 



Figure [2] shows the experimentally determined DC bias versus the root-mean-square 
driving potential, as well as modeling results for various pressures. The dashed line 
is a linear fit to the experimental data. The results from the model, calculated by 
adjusting the DC potential until the electron and ion current to the powered electrode 
balance, fall on the line, with the DC bias becoming more negative for lower pressures, 
which is in agreement with the literature O [23] , and discussed below in section [Jj 

4-2. Absorbed input power. 

The electrode power is experimentally determined by measuring the (root-mean- 
square) potential applied to the powered electrode, V rms and the current reaching 
the electrode, /. Assuming that for the given powers and pressures examined 
the plasma acts much like a simple impedance, the total power dissipation in the 
plasma is then P<n s — V rms I 24J. In the model, the total power absorbed by the 
charged species is calculated from the volume integral of the Ohmic heating term, 
P a bs — f f f (Je ■ E + J + ■ E) dV. The resulting power versus driving potential plot 
is shown in figure 03 The dotted and dashed lines are second order polynomial fits 
to the data, with coefficients of determination, i? 2 , greater than 0.99, which implies 
the fits are highly reliable. We observe that the modeled absorbed power P a i> s and 
the measured dissipated power Pdis differ by a factor of 2-4, which might indicate a 
difference in the model and actual plasma conductivity, as discussed in section [7J 




(10) 
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Figure 2. The experimentally determined DC bias (open squares) together with 
a linear fit to the data (dashed line). The circles are results from the model, 
showing good agreement with the experiment as well as a dependence on the 
pressure in agreement with the literature. 
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Figure 3. The dissipated power, P^is, measured experimentally (closed symbols, 
'E') and the total absorbed power, P a b s , computed by the model (open symbols, 
'M') and second order polynomial fits versus driving potential, for different 
neutral pressures. The fits are very reliable (R 2 > 0.99), which shows that 

Pabs,Pdis OcC(V- Vb) 2 . 
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5. Results for local parameters 

Dust clouds will be levitated in the modified GEC cell under the influence of different 
forces. We consider the electrostatic force, which depends on the dust charge and 
the electric field, and the ion drag force, which depends on the drift velovity of the 
ions, but mainly on the ion density in the discharge. The equilibrium position of dust 
particles is determined by the position where the net force acting on the dust particles 
vanishes. 

In the vertical direction, this balance is mainly determined by the electrostatic 
force and gravity (depending on the size of the particles, the ion drag might play 
an important role in the force balance). Since the vertical electric field required 
to balance gravity, E = mpg/qDi is rather large for micrometer sized melamine- 
formaldehyde particles, the particles will in general levitate near, or in the sheath 
above the powered electrode. 

In the radial direction, the electric field due to the plasma potential is generally 
relatively weak. In order to create a potential well in which the particles will be 
trapped, an additional electric field is created by the cutout in the lower electrode. 
This additional electric field causes an additional radial force pulling the particles 
inward. Usually a balance between this electric field and the inter-particle interaction 
is assumed to determine the radial equilibrium of the particles. However, it has been 
shown in micro-gravity experiments that the ion drag force is responsible for the 
formation of a dust-free void in the plasma bulk, and has thus been shown to have a 
rather large magnitude. It is therefore not unlikely that the ion drag force plays an 
important role in the radial force balance. 

The probe measurements were taken at the midplane between the two electrodes, 
which is above the equilibrium levitation height for most particle sizes, however, the 
profile behavior for varying discharge parameters will still give a good idea of the effect 
of these parameters on the forces acting on dust particles suspended in the discharge. 
These measurments can not only be used to compare the profiles obtained from our 
model with the experiment, but also to show differences between the obtained results 
and similar results reported for the original GEC reference cell. Since we were only 
able to perform Langmuir probe measurements at different radial positions, we can 
not make any statements on the change of the forces with height. Therefore, we 
mainly focus on the radial forces acting on dust particles that would be levitated in 
the discharge. 

In this section the different radial plasma profiles obtained with the Langmuir 
probe at different applied powers, and at pressures of 200 mTorr and 400 mTorr, are 
compared to similar profiles obtained with the model at the same driving potential, 
determined from figure [3l This means, for instance, that the experimental results at 
400 mTorr and 6 Watts, which are run at roughly V rms = 90 V are compared to a 
model run at the same pressure and driving potential, which is roughly 2-2.5 Watts of 
Ohmic heating. In the figures, the labels in the legends refer to the pressure and the 
power, so that 200 — 6 refers to results from a discharge at 200 mTorr and 6 Watts of 
power. 

5.1. Dust charge. 

The equilibrium charge on a spherical particle depends on the size of the particle 
and the electron temperature. The Langmuir probe determines the local electron 
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The electron temperature determines 
^ |25| . while the dust charge number 



temperature, as well as the floating potential, 
the dust charge via Z D = -\T e (eV) In f;f^ 

can be determined from the measured floating potential through Zd — Vfi(ATreoa/e). 
Because of the larger error in the electron temperature obtained from the Langmuir 
probe, we use the second approach. Figure[3]shows the dust charge number determined 
from the measured floating potential for particles with a radius of 4.45 /xm, which are 
often used in experiments in the modified GEC cell, at different radial positions above 
the lower electrode, as well as results for corresponding cases obtained with the model. 



4.0 

J> 3.5 



-i b 3.0 

N 



o 2.5 
E 



2.0 



1.5 



% 10 
& 0.5 



■ 200-1 

• 200-6 
200-1 (M) 



4.0 

J> 3.5 
o 

-i i, 3.0 
N 

£ 2.5 
E 

= 2.0 

<u 

ro 
.c 

u 

« 10 

g 0.5 



10 20 30 40 50 
Distance from center (mm) 



60 



□ 


400-1 


o 


400-6 




400-2.5 (M) 




10 20 30 40 50 
Distance from center (mm) 



60 



Figure 4. The surface charge for particles with a radius of 4.45 £tm. The results 
for 200 mTorr are shown on the left (a), the results for 400 mTorr on the right 
(b). The dashed lines represent model results for 200 mTorr at 1 Watt, and 400 
mTorr at 2.5 Watts. 



The dust charge profiles differ significantly for the low power (1 Watt) and high 
power (6 Watts) case. For low power the dust charge shows a strong dip around 20 
mm for both pressures, followed by a peak at 30 mm from the center, whereas at high 
power, the dip is insignificant and the dust charge decreases smoothly outwards. It 
is also clear that the dust charge decreases for higher pressure, consistent with our 
model results, and increases slightly with increasing power. A similar reduction of the 
dust charge with increasing pressure was also reported in [26j . The difference between 
the model and the Langmuir probe data is roughly a factor of 2, which we expect has 
to do with the ion-neutral collisions. (See the next section and a discussion in section 
0) 



5.2. Plasma densities. 

The electron and ion densities in the experiment are shown in figure [5] and figure [6l 
respectively, together with model results. The general trend for the plasma densities 
is to increase with the absorbed power, however, the response of the ions and electrons 
to a change in neutral gas pressure is different. 

We observe that the electron density increases roughly 6-fold when the dissipated 
power increases 6-fold, which might indicate that n e cx P a b s , even though there are 
only two values of the power considered here. The electron density decreases with 
neutral pressure, which is a result quite different from other results reported in the 
literature. The shape of the radial electron density profiles corresponds very well with 
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Figure 5. The electron densities measured in the discharge (symbols) and two 
corresponding model solutions (dashed/dotted lines). Note that the model values 
have been divided by arbitrary factors, showing that the absolute value of the 
electron density predicted by the model is much higher. 



the profiles obtained in the model, however, the absolute value of the electron density 
is much higher in the model. We also see that the difference between the model and 
the experiment changes with neutral pressure. This decrepancy might have to do with 
the decrepancy we observed between the absorbed and dissipated power. This will be 
further discussed in section [7] 
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Figure 6. The Ar+ densities at 200 mTorr on the left (a) and at 400 mTorr 
on the right (b) measured in the discharge (symbols) and the model solutions 
indicated by the dashed lines. The corrected ion densities are explained in section 
[7] and have to do with ion- neutral collisions. 



The ion density increases with power, but by a factor 2-3 rather than a factor of 
6. For increasing pressure, we see that the ion density decreases, but only by a small 
fraction. The ion density profile peaks in the center and falls off to the edge of the 
electrode. The model profiles are similar to the electron density profiles, having a peak 
off-axis. Furthermore, the model ion density is much higher than the experimentally 
obtained ion density. 
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A caveat is in place here, however. It has been shown in the past that the ion 
density is underestimated by ordinary Langmuir probe data in the original GEC cell 
at pressures above 50 mTorr [J]. It was shown that for 50 mTorr < P gas < 250 
mTorr, the ion density is roughly a factor of 2 higher than the value obtained from 
a Langmuir probe, whereas for P gas > 250 mTorr, it becomes a factor R pr obe/lmfp 
higher, with R pro be the radius of the probe-tip and l m fp the ion-neutral momentum 
transfer mean-free path. At 400 mTorr, this results in a factor of 5. We have added 
these corrected ion density profiles in the graphs. They are reasonably in agreement 
with the model results. 

Finally, for all powers and pressures, the ion density is higher than the electron 
density, so that the modified GEC cell seems to be strongly electropositive, which 
might indicate that a clear quasi-neutral bulk does not exist for these discharge 
settings. This is also found with the model, but the difference between the electron 
and ion density is much smaller. 



5.3. Plasma potential. 

Figure [7] shows the measured plasma potential profiles as well as profiles found from 
the model, for 200 mTorr and for 400 mTorr. We see that the plasma potential shows 
the same behavior as the plasma densities, i.e. an increase with power, but a decrease 
with pressure. The model results in this case are close to the measured profiles, even 
though the model always results in an off-axis maximum in the plasma potential, 
related to the off-axis maximum in the plasma densities obtained in the model. 
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Figure 7. The plasma potential at 200 mTorr (left, labeled a) and 400 mTorr 
(right, labeled b). 



Although the profiles are relatively flat, we still see a distinct difference between 
the low pressure (200 mTorr) and high pressure (400 mTorr) case, namely that the 
former has a positive radial derivative everywhere, while the latter has a negative radial 
derivative in the inner region of the discharge. This means there is a small electric field 
pointing inward for lower pressure, and hence a small electric force pointing outward 
for negatively charged particles, whereas the opposite is true at higher pressures. 
Therefore, the electrostatic force seems to add to the radial confinement at higher 
pressures, but acts against the confinement at lower pressures, at least in the plane 
between the electrodes. 
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6. Forces acting on the dust 

Even though we do not have any data that would allow us to show any axial (vertical) 
force dependence on the discharge parameters, we can use the observed dust charge to 
see how large the electric force would have to be to balance gravity. For MF particles 
with m — 4.5 • 1CF 6 m and p=1510 kg m~ 3 , we find that the force of gravity equals 
F g = mog = 5.65 -10~ 12 N. Using 15.000e as the minimum dust charge from figured! 
we find that the upper bound for the vertical electric field required to balance gravity 
is E zmax = m D g/Q D w 2400 V/m. This is roughly 10 times the value of tyical 
ambipolar electric fields, which shows that the particles will be suspended well below 
the midplane of the electrodes, towards the lower sheath. 

In the radial direction, we are interested in the electrostatic force and the ion 
drag force. The Debye length, Ad w A+ = 69y/T+/n+ can be determined from the 
obtained ion density. The Coulomb logarithm, A, can then be determined with the 
knowledge of the dust charge. The radial electric field can be determined from the 
gradient in the plasma potential, which, together with the ion density gradient, results 
in the ion drift velocity from equation [1] Knowing all these parameters, the ion drag 
can be determined. The electrostatic force is determined from the dust charge and 
the determined electric field. 

The resulting net force obtained from the Langmuir data in this way is shown 
in figure [5] for 200 mTorr (a) and 400 mTorr (b) . For the ion density we used the 
densities obtained by the Langmuir probe corrected for ion-neutral collisions. We see 
that for high pressure the force is directed inward near the center, whereas for lower 
pressure the force is directed outwards almost everywhere, except near the outer edge 
of the electrodes. 




10 20 30 40 50 60 10 20 30 40 50 60 



Distance from center (mm) Distance from center (mm) 

Figure 8. The net radial force at the midplane calculated from the Langmuir 
probe data for 200 mTorr on the left (a), and 400 mTorr on the right (b). 

The radial force profile obtained from the model is similar for both pressures, since 
no electric field reversal is apparent in the model results. The model force profile at 
400 mTorr and 2.5 Watts is shown on the right in figure [H Both the magnitude as well 
as the shape are in good agreement with the experimental results. This indicates that 
the corrected ion density, obtained with the Langmuir probe, as well as the measured 
plasma potential profile (and its derivative) are reliable, since the model equations for 
the forces are generally believed to be accurate. 
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7. Discussion and outlook 

In this section we discuss our results, point out the differences between the current 
results and the original GEC cell results reported in the literature, and discuss the 
possibilities for future work. We start our discussion with the results obtained for the 
global parameters of the modified GEC cell used for dusty plasma experiments. 

7.1. Natural DC bias. 

A first important conclusion for the model is that the difference between the model 
and the actual experimental geometry (i.e. the hollow grounded electrode vs. a 
solid electrode), thus the difference in area interacting with the plasma, is negligible, 
meaning that the model geometry is an excellent approximation to the real geometry. 
When we measure the slope of the fit in figure [2] and use the result in equation 
[TOl we obtain a ratio of A g /Ad — 4. The surface of the powered area is simply 
Ad — "KTrp ~ 92 cm~ 2 . The total available grounded area in the experiment, 
A g w 3500 cm~ 2 40^4^- Apparently, the grounded surface area interacting with 
the plasma [plasma facing) is much smaller than the total available grounded surface. 

Since the floating potential is directly related to T e , the dependence of T e on 
the pressure should explain the behavior of the DC bias with pressure. The global 
balance that determines the electron temperature is the balance between the creation 
of particles (through electron impact ionization, with a ionization constant Ki Z (T e ) 
(m 3 s -1 )) and the losses of particles (through the flux of particles to the walls at the 
Bohm velocity, ub — \JkT e /m + ). For a cylindrical discharge with radius R, height 
H, gas density n gaSl and electron density n e this balance can be written as [27] 

ttR 2 H — {2nR 2 a + 2irRHb) n e u Bl (11) 



with a = 0.8/ 'y / 4 + R/l m fp and b = 0.86/^3 + h/2l m f P - Here l m f P = l/n gas a m the 
mean-free path for ion-neutral collisions (including both clastic scattering and charge- 
exchange collisions, cr,„ s» 10~ 18 m 2 ). This form of the particle balance is valid for 
Pgas ~ 100 mTorr. Reordering the terms, we find 

Ka.4//)- 1 = (12) 

where d e // = 0.5RH/(aR + bH). We obtained a good approximation for n gas d e ff in 
terms of the pressure, P gas = n gas kBT gas , with R = 0.1 m and H = 0.2 m for our 
modified GEC cell: 

i / \ ^ \fk-BTgas ( Pgas \ ^ no \ 

n gas d eff {n gas ) w ) ' ( 13 ) 

whereas a numerical approximation of the ionization coefficient is given by |28j : 

K lz = 4.8 • 10- 17 (^-T e (eV) - 5.3] exp ( , ~ 4 ' 9 I . (14) 



2 



5.3 



Using a Taylor expansion around T e ps 5 eV, which is the electron temperature we 
expect in the modified GEC cell, we can find an approximation to the ionization 
coefficient up to second order, which is correct within 4% for 4 eV < T e < 7 eV: 

K lz {T e ) « 3.9 • 10" 18 + 7 • 10" 18 (T e - 5) + 3.2 • 10~ 18 (T e - 5) 2 . (15) 



Experimental and computational characterization of a modified GEC cell. 



15 



With u B {T e ) w 1550-s/TjeV) a good approximation for K iz (T e )/u B {T e ) valid for 4 
eV < T e < 7 eV is given by 

^-( T e) _ 2J . 1Q -21 (^e - 3.9) 2 (16) 



u B (T e ) ■ v?; 

Using equation (fT5|) and ([TB|) in equation (fT2")l and taking the logarithm on both 
sides we find: 



log (i^)^ log ( T " /3 ) ^ T ^i%- (17) 

It follows that the electron temperature decreases with increasing pressure, but the 
change in electron temperature for 200 mTorr < P gas < 400 mTorr is small. This 
means that the floating potential at a given driving potential becomes less negative 
for increasing pressure, but that the effect is small, which is also the result obtained 
with our model. 

The above behavior for the DC bias is very similar to the behavior in the 
original GEC reference cell, and was found both in modeling [3], as well as in many 
experimental studies, for instance in [35] • We therefore conclude that the modifications 
in the modified GEC cell studied here have no significant effect on the global behavior 
of the floating DC bias. 

7.2. Dissipated and absorbed power. 

By measuring the plasma current and the applied potential, the plasma impedance 
was determined at different pressures. The results are shown in figure [9] on the left. 
The data is perfectly fitted by an exponentially decreasing impedance. Since the 
conductivity, a, of the plasma is the reciproke of the impedance, the root-mean square 
absorbed power, Pdis = o~E RMS can be determined, if we assume a certain electric 
field. Figure [5] on the right shows the dissipated power assuming an electric field, 
Erms = 800 V m _1 (represented by the blue crosses). 

The power dissipated in the experiment at Vrms = 70 V, determined form figure 
[2] is shown by the black squares. We see that the assumption of a constant root- mean- 
square electric field of 800 V m _1 at this driving potential is acceptable. Included 
in the figure is a linear fit to the experimental data, represented by the red dotted 
line, which is reasonable, and an exponential fit, represented by the green dashed line, 
which is very good, with R? > 0.99. 

Finally, the absorbed power calculated by the model at 70 V is represented by 
the open circles. The absorbed power can be perfectly fitted by a linear dependence 
on the pressure. In our model, we use an approximation for the pressure dependence 
of the electron-neutral momentum-transfer frequency, and hence the electron mobility 
EH: 

n - ( 1000 T gas \ 

Since a = en e [i e , and /i e . m oc Pgas m our model, we must have n e oc Pg as , which 
is consistent with the faster than linear increase in electron density found with our 
model. 

Clearly, there is a difference between the actual electron-neutral momentum 
transfer frequency (and hence mobility, resistivity and conductivity) and the 
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Figure 9. A) A plot of the impedance determined in the modified GEC 
cell. B) The dissipated power at Vrms = 70 V for the experiment (black 
squares), together with a reasonable linear fit (red, dotted line) and a much 
better exponential fit (green, dashed line). The absorbed power in the model 
is represented by the circles, an excellent linear fit is represented by the black, 
dashed line. 



approximation used in the model. This can explain the difference in the electron 
densities obtained with the model and the ones measured in the modified GEC cell. 
Obviously, the determination of a better approximation for the electron mobility in 
the model will be part of future work. 

The fact that the dissipated and absorbed power are best fit by a second order 
polynomial of the driving potential, rather than just a purely quadratic dependence, 
is due to the non- vanishing electric field in the electropositive discharge present in the 
modified GEC cell. We can write the electric field as a time-averaged part and an 
oscillating part (E = Eq + E), so that we find for the time-averaged absorbed power 

P ~ oWdV ot El + 2E E + E 2 . (19) 

In a symmetric discharge, the average electric field, Eq, is very small. Since E oc 
V(t) oc Vrf s'm(ojnFt), we would simply have E ■ E ~ 0-5 1 ^E"! 2 oc 0.5V^ F . However, in 
the modified GEC cell, the time averaged space charge does not vanish, as seen in the 
results for the plasma densities in figure [5] and figure [6j and the time averaged electric 
field is considerable. 

The behavior obtained for the absorbed power is similar to the observed pressure 
dependence of the power in the original GEC cell [31]. The modifications in our 
modified GEC cell therefore do not significantly change the global behavior of the 
absorbed power strongly. We now proceed by discussing the local parameters. 

7. 3. Local parameters. 

The observed dust charge number lies in the range 5 • 10 3 < Zo < 4.2 • 10 4 , increases 
with power and decreases with pressure. The last observation can in part be attributed 
to the effect of the pressure on the electron temperature, similar to the effect this has 
on the DC bias. A measurement of the dust charge, using single dust particle-dust 
particle interactions in the originial GEC cell at low pressures 32J indicated that the 
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dust charge number for similar sized particles was in the 10 4 - 10 5 range for screened 
charged particles, whereas it was always lower than 4 • 10 3 for unscreened particles. 
The first assumption seems to be consistent with our measurements. 

In [33j the dust charge was shown to decrease with pressure at pressures similar 
to our study. This was determined from the corresponding decrease in coupling 
parameter, T, for increasing pressure. It was also observed that the coupling decreased 
for increasing power. We conclude that this must be caused by a change in the 
screening length, due to an increase in plasma density, or due to an increase in the 
dust thermal motion, since the dust charge actually increases with power. 

The shape of the electron density profile in the experiment is well represented 
by the model, but the ion density profile measured in the modified GEC cell is very 
different from the measured electron density profiles, having its maximum in the center 
and falling off towards the electrode edge, rather than having an off-axis maximum. 
This results in a large positive space charge in the center of the discharge, while 
the space charge is more negative towards the edge. The model also predicts higher 
ion densities, resulting in an electropositive space charge, but the difference between 
electron and ion density is much smaller in the model than found in the experiment. 
A similar study in the original GEC cell [3] showed that the electron density and the 
ion density both have the same off-axis maximum, which is distinctly different from 
our finding. It did conclude that the ion density was higher than the electron density, 
so that the original GEC cell was also found to be electropositive at similar settings. 

The plasma potential profiles obtained in the model are close to the measured 
profiles. The difference lies in the slope of the profiles. For higher pressures, the 
measured profiles result in electric fields pointing radially outwards, while lower 
pressures result in profiles pointing inwards. The model only gives profiles resulting 
in inward pointing electric fields, because for all pressures the net space charge has a 
clear positive maximum value towards the electrode edge. 

For dust particles immersed in the discharge, the observed change in the net 
radial force for different powers would mean an increase in radial compression for 
higher pressures. It is quite likely that the change in net force (even though the net 
forces are small) could be seen as a decrease in the inter-particle distance at higher 
pressures, due to the increased radial compression. Even though it is not directly 
obvious that the measured profiles remain the same towards the powered electrode 
below, it is likely that the radial confining electric field induced by the cutout in the 
electrode is affected by the change in the plasma radial electric field, due to a change 
in gas pressure. Therefore, the shape of a suspended dust cloud and the amount of 
dust particles that can be confined at a given height above the radial cutout might 
very well change dramatically for different neutral gas pressures. 

In short, the biggest discrepancy between the model and the experimental results 
are the plasma density profiles. It is likely that this discrepancy has the same basis as 
the discrepancy between the absorbed and the dissipated power, and has to do with the 
behaviour of the electron-neutral momentum transfer and its behavior with pressure. 
Overall, the local parameters seem to be close to the local parameters determined in 
the original GEC cell. Of special note, the change in the radial dust charge profile for 
a change in pressure, and the different shape of the ion density profile might be due to 
the change in geometry with respect to the original GEC cell. This means that despite 
the excellent agreement between the modified cell and the original cell when global 
parameters are considered, differences in the local profiles are still to be expected. 
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7.4- Future work. 

The linear dependence of the model mobilities on neutral gas pressure is correct for 
low pressures, P ga s < 10 mTorr, but a transition occurs for slightly higher pressures 
|34j . which means that an updated form for the electron-/ion-mobility is required. 
Using the obtained power dissipation data to obtain a better fit for the mobilities for 
a wider range of pressures is a possibility. 

Obtaining two- or three-dimensional profiles of the plasma parameters, especially 
above the cutout on the lower electrode, is a logical next step in order to further 
benchmark our model, as well as to better understand the force balances for particles 
suspended in the parabolic confinement potential, especially for varying pressures and 
powers. These measurements could be made with a Zyvex S-100 nanomanipulator 
currently in use on another plasma cell in the CASPER laboratory, which can be 
equipped with different types of probes and can be very accurately steered throughout 
a large volume of the discharge [35] . 

The results from the discussed fluid model can be used as input and/or boundary 
conditions for MD models of collections of dust particles, that take the interaction 
potential between the charged particles into account. In this manner, dust crystal 
properties, such as the coupling parameter and bond angle, could be computed 
not only for random input parameters, but for plasma parameters consistent with 
observations. Phase transitions observed in MD simulations and experiments could 
thus be linked to the plasma parameters through the use of the fluid model, a good 
example of which is shown in [36] . 

8. Summary 

We used a 2D fluid model to simulate a modified GEC reference cell used for dusty 
plasma experiments. We obtained the global characteristics in the form of the 
dissipated power and the natural DC bias. The model geometry is in excellent 
agreement with the modified geometry, showing that only a small fraction of the 
grounded surface area interacts with the plasma. The absorbed power versus root- 
mean-square potential behavior is well reproduced, but a difference in the behavior 
versus pressure indicates that an updated form for the charged particle mobilities is 
necessary. Such an update could be reproduced from a fit of the absorbed power for 
different pressures, which should give the plasma impedance. 

The local dust charge, plasma densities, plasma potential and radial force profiles 
were obtained from the data and the model. The measured dust charge was higher 
than the dust charge obtained with the model, by a factor of 2. The dust charge profile 
shows a distinct dip and peak near the edge of the cutout for higher pressures. The 
ion density profile is quite different from the electron density profile, with a maximum 
value in the center and falling off towards the outer edge. In the original GEC cell, the 
ion density profiles were reported to be similar in shape to the electron density profile. 
The plasma potential profile at the midplane changes such that a radially confining 
electric field is produced at higher pressures. This might very well change the confining 
electric field closer to the powered electrode as well, resulting in a significant change 
in the radial confinement and the dust crystals formed in the sheath above the cutout. 

Overall, our study shows that similar behaviour of global parameters, usually 
easily measured from externally accessible electronic signals, does not necessarily 
guarantee similarity of the local parameters, which are usually only accessible by 
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measurements in-situ. This means that different modifications to the GEC reference 
cell destroy the standardization, meaning that for every modified GEC cell the 
pressure-power behavior has to be catalogued once again. 
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